Critical domain-wall dynamics of model B 
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Abstract 

With Monte Carlo methods, we simulate the critical domain-wall dynamics of model B, taking 
Q-i the two-dimensional Ising model as an example. In the macroscopic short-time regime, a dynamic 



scaling form is revealed. Due to the existence of the quasi-random walkers, the magnetization 
shows intrinsic dependence on the lattice size L. A new exponent which governs the L-dependence 
of the magnetization is measured to be a = 0.243(8). 
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According to Hohenberg and Halperin jl|, a dynamic system with no conservation laws 
is called model A, while that with a conserved order parameter is model B. For critical 
systems, the equilibrium state of model B is in a same universality class of model A. For 
lattice models, such as the Ising model, the critical temperature of model B is also the same 
as that of model A. However, the dynamic universality class of model B is different from 
that of model A. For the two-dimensional (2D) 4 theory, one has derived that the dynamic 
exponent z = 4 — rj, with rj being the well known static exponent In other words, z is 
not an independent one, and its value is around 4, much larger than z ~ 2 of model A. 

The dynamic scaling behavior of model A has been understood around and even far from 
equilibrium [l Based on the dynamic scaling form in the macroscopic short-time regime, 
new methods for the determination of both dynamic and static critical exponents have 
been developed Recent progress in the short-time critical dynamics includes, for 

example, theoretical calculations and numerical simulations of the XY models and Josephson 

fin fin fin 

8Hllj. magnets with quenched disorder [12I414I]. ageing phenomena [15Nl7j. 



junction arrays 
and various applications and deve 
domain wall has been concerned 
at zero or low temperatures 
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opments 18M22I]. Very recently, critical relaxation of a 



241 ] , and it is also relevant for the domain- wall dynamics 



The dynamics of model B is important in various fields. For example, the phase order- 
ing dynamics of model B describes the spinodal decomposition of binary alloys and phase 

|. The dynamics of driven lattice gases also belongs to model 



separation o 
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fluids [27, 



311 ]. However, the critical dynamics of model B is less studied in the literature, 
compared to that of model A 15|] . This is mainly because its critical slowing down is more 
severe, and it is difficult to reach the equilibrium. From this view, it is instructive to explore 
the short-time critical dynamics of model B, for the critical slowing down does not disturb so 



much the simulations 



In Ref . 



291 ]. the short-time dynamic scaling form is applied 
to numerically identify the universality classes of anisotropic driven lattice gases. Since the 
dynamic system is rather complicated, it induces a controversy 32l. l33|. 

Our thought is to clarify the short-time dynamic scaling form of model B starting from 
simpler systems. For the 2D Ising model, critical relaxation with a disordered initial state has 
been simulated with the Kawasaki algorithm 3414371] . The results support that the dynamic 



exponent is z = 4 — rj = 15/4. To fully understand the short-time critical dynamics of model 
B, we should explore the dynamic effects of different initial conditions. An ordered state can 
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not be the initial state of model B. But the dynamic relaxation with a semi-ordered initial 
state is important. In fact, it describes the dynamic evolution of a domain wall. In the case 
of model A, there emerge plenty of new phenomena 
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The purpose of this paper is to investigate the critical domain-wall dynamics of model 
B with Monte Carlo methods. To be specific, we simulate the dynamic relaxation of the 
2D Ising model starting from a semi-ordered state with the Kawasaki algorithm. The semi- 
ordered initial state consists of two fully-ordered domains with opposite spin orientations. 
As time evolves, the domain wall between two domains roughens, and looks like a growing 
interface, i.e., the so-called domain interface. Such a dynamic process is inhomogenous in 
space, very different from that with a disordered initial state. 

In Monte Carlo simulations, the 2D Ising model is defined on a rectangular lattice, with 
a linear size 2L X in the x direction and L y in the y direction. Anti-periodic and periodic 
boundary conditions are adopted in the x and y directions respectively. Initially, spins are 
taken to be positive on the sublattice L x x L y at the right side and negative at the left 
side. With the semi-ordered initial state, we update the spin configuration by exchanging 
two spins in the nearest neighbor with the heat-bath algorithm at the critical temperature 
T c , up to a maximum time tu = 10 6 . Most simulations are performed with L x = L y = L. 
The total of samples for average is 30 000, 20 000, 3000 and 2000 for L = 32, 64, 128 and 
256 respectively. Statistical errors are estimated by dividing the samples into two or three 
subgroups. If the fluctuation in the time direction is comparable with or larger than the 
statistical error, it will be taken into account. For convenience, we set the x-axis such that 
the domain wall between the positive and negative spins is located at x — 0. So the x 
coordinate of a lattice site is a half-integer. 

For simplicity, we first set L x = L y = L. Due to the semi-ordered initial state, the time 
evolution of the dynamic system is inhomogeneous in the x direction. Therefore we measure 
the magnetization function of x and t, 



M(t,x,L) = - 



L 

^ &xy if) 
y=l 



(1) 



Here S xy (t) is the spin at the time t on the lattice site (x,y), and < . . . > represents the 
statistical average. Generally, the magnetization may depend on L. 

At the critical temperature, there are three spatial scales in the dynamic system, i.e., the 
non-equilibrium correlation length £(£), the coordinate x and the lattice size L. One may 



believe that £(t) grows in a universal form ~ t x l z in all spatial directions, because of the 
homogeneity of the interactions in the Hamiltonian. Therefore, standard scaling arguments 
lead to the scaling form of the magnetization 

M{t, x, L) = t^ /uz M{t 1/z /x, t 1/z /L), (2) 

where /3 and v are the static exponents, and z is the dynamic exponent. For the 2D Ising 
model with the dynamics of model B, theoretical values of the exponents are (3 = 1/8, v — 1 
and z = 15/4. On the right side of the above equation, the overall factors t~^l vz indicates the 
scaling dimension of M, and the scaling function M(s,u) describes the scale invariance of 
the dynamic system. In the macroscopic long-time regime, i.e., the regime with > L, it 
has been well known that the dynamic system exhibits a universal dynamic scaling behavior. 
In this paper, we assume that the scaling form in Eq. ([2]) holds already in the macroscopic 
short-time regime, i.e., the regime with £(£) <C L, after a microscopic time scale t m , c . 

For the dynamics of model A, the dynamic scaling behavior of the magnetization is 
relatively simple. In the short-time regime, the finite-size effect of M(t, x, L) is negligible, 
because of £(£) <C L. M(t,x,L) — > t~P/ uz F(s) with s = t 1 ^ jx. For a large s, i.e., inside the 
domain interface, F(s) exhibits a power-law behavior F(s) ~ s~^°^ v with (3o/v ~ 1. For a 
small s, i.e., outside the domain interface, F(s) ~ const [23|. In the long-time regime, the 
magnetization may depend on the lattice size L, typically dominated by an exponential law 
M(t,L) ~ exp(—t/L z ). For the dynamics of model B, the dynamic scaling behavior of the 
magnetization is rather complicated. It is anomal that the magnetization shows intrinsic 
dependence on the lattice size L even in the short-time regime. In Fig. [1] (a), the time 
evolution of the magnetization is displayed. For a fixed x, M(t, x, L) obviously varies with 
L. Although M(t, x, L) strongly depends on L, it is insensitive to the boundary condition. 
In Fig. [1] (a), numerical simulations with periodic and free boundary conditions in the x 
direction are also included. All boundary conditions lead to the same results. 

The domain-wall motion of model B is driven by interchanging positive and negative 
spins in neighbor. When a negative spin escapes from the domain interface and jumps into 
the positive domain, it becomes quasi-free and moves randomly. We call this spin 'a quasi- 
random walker 1 . The average moving distance of the quasi-random walkers in a time t is the 
order of l r (t) ~ y/i. Even in the short-time regime, i.e., the regime with ~ t l / z <^ L, 
it may occur L < l r (t) ~ for the dynamic exponent z = 15/4 > 2. Therefore the 
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quasi-random walkers easily touch the boundary of the dynamic system. This should be the 
physical origin of the L-dependence of the magnetization in the short-time regime. 

In fact, the average moving distance l r (t) of the quasi-random walkers is an additional 
spatial scale for the dynamics of model B. In the literatures, such a diffusion length scale 



l r (t) ~ t^ z wit 
initial state 



wit 
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z/2 has been detected in the dynamic relaxation with a disordered 



variable l r (t)/L. The dynamic scaling form in Eq. (j2J) holds only when the quasi-random 
walkers reach a 'homogeneous' state, i.e., l r (t) 3> L. In other words, the quasi-random 
walkers should touch the boundary, turn back to the domain interface, and then stabilize 
at a homogenous state. At x — 0.5, for example, the shortest time ts for this movement is 
about ts — 4:L 2 in our numerical simulations. 

Since M(t, x, L) relies on t through two scaling variables in Eq. ([2]), its dynamic behavior 
is complicated. Nevertheless, a dynamic scaling form is typically characterized by power-law 
behaviors. Let us now concentrate on the following features of the scaling function M(s, u), 

~ r G(u)s-M" s^oo 

M(s,u)^{ (3) 
F(s)u a s — y oo, it — y 

with s = t x l z fx and u = t l l z / L. Actually, F(s) = s~^°^ u . We use the notation F(s) for 
convenience. The exponent fi^jv characterizes the spatial behavior of the magnetization 
inside the domain interface, which appears similar to the surface exponent defined on a 



surface 
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241 ]. Thus we call it the interface exponent, a is a new exponent exclusive 



for the dynamics of model B, for describing the finite-size dependence of the magnetization. 

In Fig. [1] (b), the magnetization is plotted versus x. Inside the domain interface, i.e., for 
x < ^(t) or a large s, M(t,x,L) exhibits a power-law behavior. The slope of the curves is 
0.991(8), very close to 1. This suggests (3 /u = 1, and indicates that M(t, x, L) is an analytic 
function of x, the same as that of model A 23J]. In Fig. [1] (a), we detect a power- law decay 
of the magnetization inside the domain interface in the short-time regime, i.e., for large s 
and small u. The slope is estimated to be 0.362 from the curve of L = 128. For L = 64, it is 
reaching the long-time regime for t > 10 5 , and a deviation from the power law is observed. 
This power-law behavior indicates (/3/i/+/3o/i/+a)/z ~ 0.362, and yields a ~ 0.242 with the 
input /3/u = 1/8 and z = 15/4. Since a is positive, the magnetization in Eq. (j3]) increases 
with L. Therefore we call this L-dependence of the magnetization intrinsic. 
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To illustrate the scaling form comprehensively and precisely, we now perform scaling plots 
with Eq. ([2]). In Fig. [2] (a), the scaling function M(s, u) with a fixed u is displayed, and data 
collapse is observed for different L. In the large-s regime, i.e., inside the domain interface, 
M(s, u) exhibits a power- law behavior s~^/ v with p /u = 0.99(2). This agrees with the 



measurement in Fig. d](b). 
fitted to the error function 



or the dynamics of model A, the scaling function M(s) can be 
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241 ] . For the dynamics of model B, M(s,u) with a fixed u 



also coincides with the error function in the large-s regime. However, it does not approach 
a constant in the small-s regime, due to the existence of the quasi-random walkers. 

According to Eq. 03]), M(s,u) ~ F(s)u~ a for large s and small u. In the inset of Fig. 
0(a), the scaling plot for M(t, x,u)t /3 ^ uz u a ~ F(s) is performed. Data collapse is observed 
inside the domain interface. With the input (3/v = 1/8 and z = 15/4, we extract the 
exponents a = 0.243(8) and p /v = 0.990(8), consistent with 0.242 and 0.991(8) estimated 
from Fig. [TJ (a) and (b). 

To explicitly reveal the L-dependence of the magnetization in different time regimes, we 
fix s and plot the scaling function M(s, u) versus u. Inside the domain interface, i.e, for a 
large s, M(s,u) ~ G(m)s _1 . In Fig. |2] (b), the scaling function M(t,x,u)t l3 ^ uz s ~ G{u) is 
displayed at x = 0.5. In the time regime t > t$ ~ 4L 2 , data collapse is observed. For a small 
u, there emerges a power-law decay G(u) ~ u~ a , and o is estimated to 0.25, in agreement 
with 0.243(8) extracted from Fig. [2] (a). For a large u, G(u) is governed by an exponential 
law. For comparison, a scaling function G(u) of model A is schematically shown in Fig. [2] 
(b). In the time regime t < ts, the scaling form is violated due to the quasi-random walkers, 
and data of different L do not collapse. To fully understand the dynamic behavior in this 
time regime, we need to include the additional scaling variable l r (t)/L in the dynamic scaling 
form. A detailed description of this kind will be presented elsewhere. 

In the inset of Fig. [2] (b), the scaling function of model B is plotted for a small s, 
i.e., outside the domain interface. Qualitatively, G(u) is similar to that inside the domain 
interface, but the effective a decreases with s. 

Finally, we should point out that the finite-size dependence of the magnetization in the 
domain- wall motion of model B is anisotropic in spatial directions. In fact, the intrinsic 
L-dependence of the scaling function M(s, u) in Eq. ([3]) does refer only to L x , not L y . Since 
the domain wall is oriented in the y direction, the random walkers does not induce anomal 
dependence of the magnetization on L y . In the short-time regime, i.e., the regime with 
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£(t) <C L y , the magnetization is independent of L y . In Monte Carlo simulations, we have 
confirmed this by fixing L x and changing L y from 64 to 128 and 256. 

In summary, we have simulated the critical domain-wall dynamics of model B with Monte 
Carlo methods, taking the 2D Ising model as an example. A short-time dynamic scaling 
form is revealed, and the scaling function is carefully computed. Due to the existence of 
the quasi-random walkers, the magnetization intrinsically depends on the lattice size even 
in the short-time regime. This is very different from the case of model A. The new exponent 
a which governs the L-dependence of the magnetization is measured to be 0.243(8). The 
interface exponent %jv takes the value 0.99(1), close to 1, the same as that of model A. It 
is a challenge to derive the scaling form in Eq. ($2§ and extract the exponent a in Eq. (J3]) 
with renormalization group methods. 
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(a) (b) 
FIG. 1: (a) The time evolution of the magnetization is plotted for different x, and L = 64, 128 and 
256 (from below) on a double-log scale. An anti-periodic boundary condition (APBC) is adopted 
in the x direction. For comparison, results of x = 0.5 and L = 128 with periodic and free boundary 
conditions (PBC and FBC) are also included. The dashed line shows a power-law fit. (b) The 
magnetization obtained with L=256 is plotted versus x for different t on a double-log scale. Errors 
of the data points are less or about one percent. The dashed line shows a power-law fit. 
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FIG. 2: (a) The scaling function M(t,x, L)t p / vz = M(s,u) with s = t 1 / z /x and u = t 1 / z /L is 
plotted versus 1/s for fixed u = 0.53 and 0.621. Data collapse is observed for different L. Solid 
lines represent the error function f(y) ~ J y exp(—x 2 )dx, and the dashed line shows a power-law fit. 
In the inset, a power-law behavior M(s, u) ~ F(s)u~ a is assumed in the short-time regime, and 
data of different u and L inside the domain interface collapse onto a single curve Mt^' vz u a ~ F(s). 
(b) A power-law behavior M(s, u) ~ is assumed inside the domain interface, and data of 

different L collapse onto the master curve Mt^l vz s ~ G(u). The solid line indicates the power-law 
behavior in the short-time regime. Departure of the data from the master curve in the time regime 
t <ts ~ 4L 2 is also displayed. The dashed line schematically shows the scaling function of model 
A. In the inset, the scaling function M(s, u) outside the domain interface is plotted versus u. 
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